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Entanglement transfer through the turbulent atmosphere 
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The propagation of polarization-entangled states of light through fluctuating loss channels in the 
turbulent atmosphere is studied, including the situation of strong losses. We consider violations 
of Bell inequalities by light, emitted by a parametric down-conversion source, after transmission 
through the turbulent atmosphere. It is shown by analytical calculations, that in the presence 
of background radiation and dark counts, fluctuating loss channels may preserve entanglement 
properties of light even better than standard loss channels, when postselected measurements are 
applied. 
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I. INTRODUCTION 

The distribution of quantum light through the turbu- 
lent atmosphere has attracted a great deal of attention 
since recent experiments [iHll have demonstrated the fea- 
sibility of quantum key distribution in free-space chan- 
nels. In this connection the question arises as to whether 
nonclassical properties of light can be preserved during 
its propagation in fluctuating media. Particularly, spe- 
cial interest is applied to the transfer of entanglement [5] 
through the turbulent atmosphere since this problem has 
important perspectives in quantum communications. 

The theory of classical light distributed through the 
atmosphere was established many years ago [6] . Phenom- 
ena such as beam wander, scintillations, beam spreading, 
and spatial coherence degradation have been explained 
in the framework of this description. It is based on Kol- 
mogorov's theory of turbulence, which is successfully ap- 
plied to a description of wave propagation in random me- 
dia. 

The theory of nonclassical phenomena, such as entan- 
glement, for light propagation in random media is less 
developed. First, one should mention the approach pro- 
posed in Ref. 7] for the description of the photocounting 
statistics of quantum light propagating through the tur- 
bulent atmosphere. The idea of this approach is based on 
the introduction of a stochastic intensity modulation for 
light propagating through random media. Another ap- 
proach describing nonclassical properties of light is based 
on the photon wave function [8|. It enables one to de- 
scribe a special class of entangled states of light in the tur- 
bulent atmosphere ^ . The concept of the single-photon 
wave function can be formulated in a more sophisticated 
way in terms of the single-photon Wigner function and 
then generalized to the description of arbitrary quantum 
states. This technique has been developed in 10] and 
applied to some interesting examples of quantum and 
classical light. 
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Recently, a theoretical model for the light distributed 
through the turbulent atmosphere and processed by ho- 
modyne detection has been proposed The idea is 

based on the description of random media in terms of 
a fluctuating loss channel. It has been shown that the 
turbulent atmosphere introduces, compared with stan- 
dard loss channels, additional noise into quantum states 
of light. The presence of such noise has also been dis- 
cussed recently in Ref. (T^]. Moreover, as has been shown 
in Ref. jllj , the nonclassicality of bright light fields (espe- 
cially fields with a large coherent amplitude) is more frag- 
ile against turbulence than the nonclassicality of weak 
light. In this connection special interest should be paid to 
entangled photon pairs and single photons, whose fields 
are weak enough to preserve their nonclassical properties. 

In this paper we give a systematic theoretical descrip- 
tion of polarization-entangled states of light distributed 
through the fluctuating loss channels in the atmosphere 
and processed by polarization analyzers. For such light 
we check violations of Bell inequalities [13, [l3| • Besides 
turbulence our consideration includes imperfect detection 
and noise counts originating from internal dark counts 
and background radiation. Moreover, we deal with a re- 
alistic parametric down-conversion (PDC) source of the 
radiation P, null. We apply the results of our consid- 
eration to analyzing a recently reported experiment [1], 
where light has been distributed over a 144-km atmo- 
spheric channel between two Canary Islands in a config- 
uration with both receivers being co-located at the same 
place. We show that with fluctuating loss channels one 
may observe larger values of the Bell parameter com- 
pared with corresponding nonfluctuating loss channels. 
This unexpected effect is caused by those (postselected) 
random events which are related to small losses. 

The paper is organized as follows. In Sec. |lT]we sum- 
marize the needed basic principles of Bell-type experi- 
ments. A systematic description of fluctuating loss chan- 
nels for the case of such experiments is given in Sec. IIIII 
In Sec. IIVI we consider the violation of Bell inequalities 
by light after propagation through the turbulent atmo- 
sphere, when the source generates perfect Bell states. 
The situation for a realistic PDC source is analyzed in 
Sec. El In Sec. |VT]we describe a procedure for measur- 
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ing the needed turbulence parameters. A summary and 
some concluding remarks are given in Sec. I VIII 



II. BELL-TYPE EXPERIMENTS 

Let us start with the consideration of a typical exper- 
imental setup (see Fig. [1} . The entangled photon pairs 
are emitted by the source S and then transferred through 
the turbulent atmosphere to the receiver stations A and 
B. In principle, different configurations of the considered 
experiment are possible: The source can be placed sepa- 
rately from both receivers (e.g., on a low-orbit satellite), 
it can be placed on the site of one of the receivers 
or both receivers can be situated at the same place 
(a configuration convenient for testing the feasibility of 
entanglement transfer) . At the receiver stations the light 
is collected by a telescope (or other device) and subse- 
quently directed to the polarization analyzers. Each of 
these analyzers consists of a half-wave plate (which turns 
the polarization direction by the angles and 0b at 
the receiver stations A and B, respectively), a polarizing 
beam-splitter, and two detectors. 
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FIG. 1: A typical experimental setup for checking the vio- 
lation of Bell inequalities for light transmitting through tur- 
bulent atmosphere. The source S produces entangled photon 
pairs. The polarization analyzers at A and B sites consist of 
half-wave plates, HWP, polarizer beam-splitters, PBS, and 
two pairs of detectors; -Dt^ for the transmitted signal and 
Drjj^ for the reflected signal at the A site; -Dtb for the trans- 
mitted signal and -Dr^ for the reflected signal at the B site. 

According to photodetection theory (T^, [13] , the prob- 
ability of registering one count at the detector ia = 
{Ta,Ra} by receiver A, one count at the detector ib — 
{Tb,Rb} by receiver B, and no counts at the other de- 
tectors is given by 



i^.,3(^A,^B) = T>(n|rfi«ii(°)ii£) 



(1) 



^A 7^ JA, 7^ JB, where g is the density operator of the 
light field at the inputs of both polarization analyzers. 
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(n) 
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('7iA(B)"'«A(E 



!±<£lZ_e"'''A(B)"'A(B) -^'A(B) 



(2) 

is the positive operator-valued measure for the detec- 
tor iA(B) ,18], '7iA(B) and iViA(B) are the efficiency and 
the mean value of noise counts (originating from inter- 
nal dark counts and background radiation) , respectively. 



for the detector iA(B), and : : means normal ordering. 
The photon number operator at the input of the detec- 
tor iA(B) can be written in terms of the corresponding 
annihilation and creation operators. 



«A(B) *A(E 



(3) 



These operators can be expressed via the operators of the 
horizontal and vertical modes, Sh. and a^j , respec- 

A(B) A(B) 

tively, at the input of the polarization analyzers as 

«^A(B) = «Ha(b) cos 6a(b) + avA(3) sin ^a(b) , (4) 



^Ha(b) sin6'A(B) 



"A(B) 

"va(b) <=0S ^A(B) ■ (5) 



In the case of using on/off detectors, which cannot distin- 
guish between different photon numbers, Eq. ([T]) should 
be rewritten as 



+ 00 

iAaB(^A,eB)= E Tr(n[: 



)nWn(°)n(o)^). (6) 



The correlation coefficient, which appears in the Bell 
theory [13], is given by 



where 



-Psamc (^A, Ob) — -Pdiffcrcnt {d A, Ob) 
^B) = — 77\ ^ N , „ 77{ [I) 



?A,»B 
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ff'crcnt ' 
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Ps. 
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Ta,Tb 
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0B)+PRA.nu{SA,0B) (8) 



is the probability of getting clicks on both detectors in the 
transmission channels or both detectors in the refiection 
channels, and 

^different (^'a, ^b) = -Pta.Rb (^A, ^b) + ^Ra.Tb (^A, ^'b) 

(9) 

is the probability to get clicks on the detectors in the 
transmission channel at one site and the reflection chan- 
nel at another site. According to the Clauser-Horne- 
Shimony-Holt (CHSH) Bell-type inequality [ll] for any 
local-realistic theory the maximal value of the parameter 



) fl(2) 



(10) 



should be 



B<2. 



(11) 



Quantum light may violate this inequality. In this case 
the maximal violating by the value 2\/2 is reached for 
the Bell state 



1S) = ^(11)ha 10)va |0)hb |1)vb 
+e'^-lO)Hjl)vJl)HB |0)vb) 

= 4(|H)AlV)B + e'nV)A]H)B), 



(12) 



for ^ = ^ and {e^^^^B^^A^^B^) = (O, f ' 
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III. FLUCTUATING LOSS CHANNELS 

The next problem is to derive the density operator g for 
the hght transmitted through the turbulent atmosphere. 
Similar to Ref. [ll|, where homodyne detection of such 
light has been studied, we will consider the atmosphere 
as an attenuating system. However, an important fea- 
ture of the present case is that the nonmonochromatic 
output mode is not specified by the local oscillator; it 
can take an arbitrary form depending on the distribu- 
tion of the refraction index in space. Generally speak- 
ing, it is impossible to describe such light in terms of 
a monochromatic or nonmonochromatic four-mode den- 
sity operator. However, it can be represented in terms of 
four nonmonochromatic modes with fluctuating shapes. 
Of course, such a representation is rather formal since it 
is difficult to provide phase-sensitive measurements for 
such a mode. However, it appears to be useful for evalu- 
ating observable quantities in the considered experiment. 

Let a"' and a'" be the field annihilation oper- 

"A(B) A(B) 

ators of the horizontal and vertical modes, respectively, 
generated by the source in the direction of the A{B) re- 
ceiver station. The operator input-output relations for 
the attenuating system can be written as 



Ha(B)«Ha(B, +^HVA(B)aVA(B, + 



^Ha(b)Cha(b) ' 

(13) 



«Va(b) ='^Va(B)«Va(b) +^VHA(B)aHA(B, + ^Va(b) 



A(E 



(14) 



where Th-.o, and Tv.,^. are the transmission coefficients 
for the horizontal and vertical modes, respectively, in 
the direction of the A(B) receiver stations; cl? and 

'^va(b) ^""^ operators describing the losses related to 
absorption and scattering with the absorption and reflec- 
tion coefficients -Rha(b) ^^^^ ^Va(b) ' respectively. Since 
the depolarization effect of the atmosphere is extremely 
small [6|], we may set in the following ThVaib) ~ 
^vHa(b) ~ 0' and Tha(b) ~ TVajb, = Ta(b)- Hence, 
the operator input-output relations ([T5| and can 
be simply written as 



«Ha(b, - Ta(B)<j,^^^ + i?HA(B)CHA(B)' (15) 

aVA(B, = ^A(B)<,3, + i?VA(B,cl"A(B,- (l^) 

The transmission coefficients Ta(b) a-nd the absorp- 

(B) 



tion and reffection coefficients i?H(v)A(B^ satisfy the con- 



stramts 



T, 



A(B)| 



R 



H(V)a(e 



- 1, 



(17) 



These relations can be converted into the correspond- 
ing relations between the density operator g-m of the light 
generated by the source and the density operator g-j, 
[where T = (Ta, Tb)] of the light transmitted through the 
loss channels under the assumption that the absorption 



and scattering system is in the vacuum state. For exam- 
ple, the relation between the Glauber-Sudarshan P func- 
tion [l^ of the attenuated light, P-r (oha , c^va , Q^hb i Q^vb ) i 
and the P function of the light generated by the source, 
-Pin (q^ha 1 "va 1 "hb 1 "vb ) 1 is given by [i3 



1 p ( ^ ^ ^ ^ 

ITaI'ITbI^ "\Ta ' Ta' Tb' Tb 



(18) 



Similarly, the corresponding characteristic function of the 
attenuated light, $t (/3ha i /^va - /^hb : -^vb ) : is related to 
the characteristic function of the light generated by the 



source. 



$in (/3ha,/3va,^Hb,^v 



as 



<i>T(/?HA,/3vA,/3HB,/3vB) (W) 

= <i>in (T1/3ha , T1/?VA , r^/3HB , T]^/3vb ) • 

The corresponding expressions can also be written for 
s-parametrized phase-space distribution, for normal- 
ordered moments (see, e.g., Ref. [lT|), for the density 
operator in the Fock-state representation [l^, and, in 
principle for an arbitrary representation of the density 
operator. 

The main difference of the ffuctuating loss channels 
from the standard loss channels is that the transmission 
coefficients Ta and Tb are random variables. This means 
that the density operator of the field at the input of the 
polarization analyzers, g, should be obtained by averag- 
ing the density operator with the probability distri- 
bution of the transmission coefficient (PDTC), V (T), 



9 - 



JdTV(T) g^ 



(20) 



where dT = cPTacPTb, and the integration domain is 
defined by the conditions 



D= {|Ta|' < 1, |Tb|' < l} 



(21) 



Such an approach is closely related to the idea of a ran- 
dom intensity modulation in the description of the photo- 
counting distribution for light transmitted through tur- 
bulent media [7,]. An important difference is that the 
factor of the intensity modulation was allowed to attain 
values in the range of [0, -|-cx)). In our approach this do- 
main is restricted to [0, 1], which preserves the required 
positivity of the density operator in Eq. ([^0]) . 



IV. BELL STATES 

We start with the consideration of the idealized situa- 
tion, when the source generates the Bell state . This 
will help us to better understand the nature of contribu- 
tions from noise, atmospheric turbulence, and different 
experimental imperfections. First, we consider the den- 
sity operator (and its entanglement properties) of such 



4 



a state after passing the light through fluctuating loss 
channels. Next, we include in the consideration back- 
ground radiation, dark counts, and detection losses - for 
such conditions we analyze violations of Bell inequalities. 



A. Output density operator 

The density operator of the Bell state (IT^ is simply 
written as ftn = Qe = \B) {B\. Utilizing Eqs. (fT8]) - (pO)) 
one gets that the corresponding density operator for the 
light at the input ports of the polarization analyzers is 
written as a convex combination, 

(22) 

of the following states: vacuum state Qq with the proba- 
bility 



Po = ((i-|Ta|') (i-ITbI 



(23) 



single-photon states ^ha i ^va j j ^'^'^ 9vb corre- 
sponding modes with the probabilities 



Pha = Pva = ^ (|Ta|^ (l - \Tb\^ 



Phb=Pv3 = 2(|TB^1-|TA|■ 



(24) 



(25) 



and the Bell state [cf. Eq. ([T^ ] with the probability 



In the absence of background radiation and dark 
counts, only the contribution of the density operator 
in Eq. ([22l) is postselected in the measurements. Hence 
the entanglement properties of the light for such an ex- 
periment are not destroyed by the atmospheric turbu- 
lence, at least in the absence of dark counts, background 
noise and other experimental imperfections. The result, 
nevertheless, may be significantly changed when such ef- 
fects play a major role. 



B. Imperfect photodetection 

In the presence of background radiation and dark 
counts, the postselection procedure no longer results in 
the perfect separation of the density matrix g'^ from the 
combination (|22|) . Indeed, in this case it is possible that 
one click at the receiver station A{B) originated from the 
signal is combined with a click originated from noise at 
the receiver station B{A). Also it is possible that noise 
contributes into clicks at both receiver stations. 

This situation should be analyzed by substitution of 
the density operator ((22l) in Eqs. ^ and ^ and then the 
result should be used in Eq. ([7]). The general analytical 
solution of this problem is given in Appendix]^ Here we 
consider an important special case of equal detectors and 
equal background noise at all detectors, that is, »7iA(B) ^ 
77c and Ni^^g^ = Nnc for aU za(b) in Eq. ©. 

After some algebra, the correlation coefficient ^ ap- 
pears to be equal to 



-(iTAriTBr 



(26) 



The brackets (. . .) mean averaging with the PDTC, 
V{T). 



E{eA,9B) (27) 
= 5 [— cos 26'a cos 26'b + cos (p sin 26'a sin 26'b] • 

For photon-number-resolving detectors we get 



S 



Pb [Vc + (1 - Vc) Nn 



Ps [77c + (1 - 7;c) iVnc] + 2p,Nnc [Vc + 2A^„c (1 - 77c)] + ApoN?,, 



and 



S 



P.e^^- [77c + (1 - 77c) (1 



PbS 



2iV„, 



[Vc + (1 - Vc) (1 - e- 



+ 2p, [e^nc _ 1] [^^eiVnc + 2 (e^f^c - 1) (1 - ,;c)] + 4po [e^- - 1]^ 



(28) 



(29) 



for on/off detectors. Herein, pi — Pha +Pva +-Phb +Pvb is the probability of the appearance of a single-photon state. 



The visibility V = E{^,^) (see [!1,[3]) can be simply 
obtained from Eq. ([27)) as 

V {ip) = S cos ip. (30) 

The maximal and minimal values for this parameter are 
V± = ztS. In the case of different detectors and/or 



different background radiations at these detectors, see 
Appendix 1X1 we obtain V+ ^ —V^. The parameter S 
strictly depends on both noise counts and turbulence 
properties of the atmosphere. Hence, in contrast to 
the case of perfect photodetection, the maximal value 
of the Bell parameter [Eq. (jlOl) ] depends on the atmo- 
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spheric turbulence. In Fig. [2] we consider this parameter 
as a function of the parameter cos ip, which appears in 
Eq. dSll). 




FIG. 2: (Color online) The maximal value of the Bell param- 
eter B vs. cos If, for different values of the parameter S: (a) 
1, (b) 0.9, (c) 0.8. 

The obtained analytical results can be used for an- 
alyzing a realistic experimental situation, such as that 
considered in Ref. In this case both receiver stations 
had been positioned at the same place. Consequently, 
the light pulses at receivers A and B have been propa- 
gating along the same path through the fluctuating at- 
mosphere. They are separated by a small time interval, 
which is much less than the characteristic time of the 
atmospheric fluctuations. Hence, in this case one can 
simply set 



\TaY 



(31) 



where ?7atm is the fluctuating atmospheric efficiency. We 
also suppose that this efficiency is approximately log- 
normally distributed (Tl| . 



V (??atm) 



1 



27rfT?7a 



■ exp 



1 ^Injyat 
2 



(T 



, (32) 



where 9 — — (InT^atm) characterizes the mean atmo- 
spheric losses and a (the variance of 6* = — luTyatm) char- 
acterizes the atmosphere turbulence. It is worth noting 
that this distribution can be applied only for a <^ 9. 
Realistic parameters for the detection efficiency and the 
mean atmospheric losses can be extracted from Ref. [l[: 
i]c = 0.25 and 9 — 7.7 (which corresponds to 32-dB 
single-photon atmosphere attenuation and 1.5-dB losses 
at the receiver device |21i]). 

An important conclusion is that the visibility in the 
presence of loss- fluctuating channels attains higher values 
compared with similar standard-loss channels (or with 
slightly fluctuating loss channels) (see Fig. [3]). This fact 
can be explained by contributions of random events with 



9 = — In r^atm (losses) less than the average value 9. In 
the presence of postselected measurements, this plays a 
significant role. 
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FIG. 3: (Color online) The visibility V+ vs. the mean value 
of noise counts A'nc for different values of the turbulence pa- 
rameter a: (a) 0.1, (b) 1, (c) 2. The detection efficiency is 
77c = 0.25; the mean atmospheric losses are 9 — 7.7. The 
result is equal for both photon- number-resolving and on/off 
detectors. 



PARAMETRIC DOWN-CONVERSION 
SOURCE 



The realistic sources used in experiments generate ra- 
diation states, which are more complicated than the Bell 
state mi). For example, in Refs. P, [2| one uses a PDC 
source for the generation of entangled photon pairs. For 
such a source, the contribution from the multiphoton pair 
emission is essential. In the presence of losses or/and in 
the case of using on/off detectors, all these photons con- 
tribute to the measured Bell parameter. Moreover, the 
result appears to be very sensitive to the background 
noise and dark counts. 

The state generated by the PDC source is given by 

(cf. El) 



+ CSO 



|PDC) = (coshx)"^X!^/^^*^^l^"^l*")' (33) 

n=0 

where x is the squeezing parameter, and 



(34) 



1 " 

= E 



e**^™ \n - to)h^ |to)v^ \m)^^ \n - m) 



771 — 



For n = 1, |$„) is the BeU state, cf. Eq. The 
state is Gaussian - its characteristic function of the 
Glauber-Sudarshan P function is written as 
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xexp 



tanh^ X |/3va r + ^^^^^ X |^hb |^ - tanhx (/3va/3hb + /^v^^hb) 
1 — tanh X 

tanh^ X |/?ha r + tanh' X |/3vb T " ^anhx (e^'^/^H^/^VB + 

1 — tanh^ X 



(35) 



We apply Eq. p5l) for the quantum-state input-output 
relations (|19l) and (l20t . In the resulting characteristic 
function, the variables (/?haj/5va;/3hb,/3vb) ^i'^ trans- 
formed to the variables (/3tai/3raj/3tb)/3rb) using the 
input -output relations for the polarization analyzers @ 
and ([5]) , which should be rewritten in terms of the argu- 
ments of the characteristic function. Finally, we substi- 
tute the resulting expression in Eqs. (O and (H]), which 
should be also rewritten in terms of the characteristic 
function. The obtained result is used in Eqs. ([7]) and 
for analyzing violations of Bell inequalities. For the 



Qa = a^=V (1 - ?7) (1 - tanh^ x) tanh^ x (41) 
X {772 tanh^ X - [1 + (?? - 1) tanh^ x] ^} , 



Cta.Tb = Cra,Rb=V^ tanh-^ X [l - tanh"^ x] ' 
X (1 — ry)^ tanh^ x ^ sin^ 



(42) 



probability Pi^ 



9a, 6*3) we get 



PiA,iB (^'a, 6'b) = »7ta%a%b%b (1 - tanh^ x) (36) Cta,Rb = C'ra,Tb='7^ tanh^ x [l - tanh^ x]' 



(43) 



X exp (-N„c 



CiA^'lB 



iV,;, 



^2 



CiA 

^0 



(^2 



1 



for photon-number-resolving detectors and 



Pi, 



3) = ^A^A^I 



X exp (-N„c) 



1 



1^1 - tanh^ x) (37) 



Cq + CiA + CiB + Ci 



Co + a, 



Co + Cw 



for on/off detectors. Here 



Ta 



lA.lB 

1 

c'o 



(38) 



is the total number of noise counts by the four detectors. 
The analytical form for the coefficients Co, Qa, Qb, and 
CiA.iB is given in Appendix IB] 

Consider again the case when both receiver stations 
are situated in the same place (see [l|). For simplicity 
we suppose all detectors to be equal and all transmis- 
sion coefficients to be strongly correlated [i.e. satisfying 
condition (|3ip ]. As follows from the equations in Ap- 
pendix |B1 in this case one can combine atmospheric and 
detection losses in one efficiency. 



V = ?yatm'7c- 



(39) 



We suppose this efficiency to be log-normally distributed, 
similar to Eq. (15^ . In this case, the coefficients Co, Qa, 



CiB, and CiA,iB have the form 



Co = {r;2 tanh^ x - [l + (r; - 1) tanh^ x] 



(40) 



(1 — 77)^ tanh^ X ^ cos^ 



The realistic values of some parameters can be extracted 
from Ref. [l|: ^ = 9.1 (32 dB of single-photon atmosphere 
attenuation, 1.5 dB of losses at the receiver device [llj, 
and detection losses with rjc = 0.25) and A^nc = 0.5 x 10~^ 
(200 counts/s from the dark counts and 200 counts/s 
from the background radiation in the time window of 
1.25 ns). 

In the case of photon-number-resolving detectors, the 
Bell parameter attains maximal values in the expected 
range (less than 2^/2 but larger than 2) for small val- 
ues of the squeezing parameter x (see Fig. |4]) . The pro- 
nounced minimum near tanhx = is caused by con- 
tributions from classical background radiation and dark 
counts. Similar to the case of Bell states, the fluctuating- 
loss channels preserve entanglement properties of light 
better than similar channels with standard (or slightly 
fluctuating) losses. This unexpected property is a conse- 
quence of contributions from random events with < 9 
in the postselected measurements, where 9 ~ — In 77. 

The situation becomes different in the case of using 
on/off detectors (see Fig. [5]). In the case of small fluctua- 
tions (or no fluctuations) of the transmission coefficients, 
the maximal value of the Bell parameter may attain val- 
ues near zero for the whole range of the parameter x- 
This is caused by a strong sensitivity of this type of mea- 
surement with respect to noise counts. When the tur- 
bulence parameter a becomes rather large, contributions 
from events with 9 < 9 in the postselected measurements 
are dominant. This explains the nonclassical values of 
the Bell parameter for strong turbulence. 
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FIG. 4: (Color online) The maximal value of the Bell pa- 
rameter B, obtained for photon-number-resolving detectors, 
vs parameter tanhx for the mean value of noise counts 
A'^nc ~ 0.5 X 10"", the mean total losses Q = 9.1, and the 
turbulence parameter a: (a) 0.1, (b) 2, (c) 3. 
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FIG. 5: (Color online) The maximal value of the Bell parame- 
ter B, obtained for on/off detectors, vs the parameter tanhx, 
for the mean value of noise counts A'^nc = 0.5 x 10~®, the mean 
total losses Q = 9.1, and the turbulence parameter a: (a) 0.1, 
(b) 2, (c) 3. 



VI. MEASUREMENT OF TURBULENCE 
PARAMETERS 

A complete theory of quantum light in the turbulent 
atmosphere should include a model for the explicit form 
of the PDTC. Of course, the log-normal distribution ([5^ 
cannot be applied in general. Alternatively, one can con- 
sider a procedure that enables one to reconstruct the 
PDTC from independent measurements. Such a pro- 
cedure has been proposed in Ref. (jLlJ for the case of 
homodyne detection of the light propagating through 
the atmosphere. In the considered case the situation 
is rather different: Since the output nonmonochromatic 
mode has a fluctuating form, it is difficult to measure 
its phase properties. On the other hand, the phase in- 



formation is not necessary in the considered experiments 
- the corresponding equations (see for example Appen- 
dices and |B|) include only mean values of the functions 
of ITaI = "Qk and |Tb| = rj^. Hence, one can consider 
the reconstruction of the PDTC, 'P(?7a,?7b), or alterna- 
tively, its statistical moments {rj^ij^)- 

Let the input light be prepared in coherent states, |q!a) 
in the direction of the receiver A and la^) in the direction 
of the receiver B. For the sake of simplicity we consider 
a realistic situation when the resulting signal at the re- 
ceivers is much stronger than the background radiation 
and the dark counts. In this case the photocounting dis- 
tribution at the receivers is given by 

1 

PuA^na ^ / drjAdriBV {r]j,,r]s) (44) 









"A 








VcVb I^bI^ 



where tjc is the detection efhciency. In principle, the 
methods of ill-posed problems (see, e.g., [22]) enable one 
to invert Eq. (gll) and to get the PDTC, P (77A, ?7b), from 
the measured photocounting statistics. 

An alternative way to resolve this problem is the recon- 
struction of statistical moments of the PDTC. A mean 
value of a function of ?7a and ?7b can be presented as a 
series of such moments. These moments can be obtained 
from the corresponding moments of photocounts, 

E «Pn.,n^- (45) 
tia.wb 

For example, for the first moments, one gets 

nX" = IckaI^ (?7a> , (46) 

|aB|'(77B). (47) 

The first statistical moments of the PDTC, {ija} and 

are easily obtained from these equations. Similarly, the 

second moments of photocounts are given by 

uaub = I^aI^ \as\^ {VaVb) , (48) 
^=\a.f{vl) + \a.f{v.), (49) 

^=\a^\''{rjl) + \a^\'{rj^). (50) 

Combining these equations with Eqs. (|46p and (|T7|) . one 
gets the second statistical moments of the PDTC from 
UAnB , n\, rig , nX, and Tzb ■ In the same way one can get 
higher moments of the PDTC by considering the corre- 
sponding moments of photocounts. 
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VII. SUMMARY AND CONCLUSIONS 

In recent experiments violations of Bell inequalities 
have been studied, by using entangled light after trans- 
mission through the turbulent atmosphere. Here we have 
presented analytical investigations of the transfer of en- 
tanglement through channels with fluctuating losses. The 
effect of the atmospheric turbulence has been modeled 
by the statistical properties of the complex transmission 
coefficient. We show how the corresponding statistical 
characteristics of such channels can be experimentally 
determined. Of course, the turbulence itself usually di- 
minishes nonclassical properties of light. In the consid- 
ered case only a small depolarization caused by the at- 
mospheric turbiilcncc might slightly destroy the entan- 
glement properties of light. However, since this effect is 
very small it has been neglected. 

For a realistic analysis of experiments, one must take 
into account that the measurements are usually per- 
formed in the presence of background radiation and dark 
counts. This may also lead, in the presence of strong 
losses, to a decrease of the measured value of the Bell 
parameter. In this context it is an important result of 
our treatment that fluctuating loss channels may pre- 
serve entanglement properties of light even better than 
the corresponding standard loss channels. The reason 
for this fact as follows: Some of the recorded events are 
related to small values of the fluctuating atmospheric 
transmission. These events are typically caused by back- 

I 



ground radiation and dark counts. On the other hand, for 
strong atmospheric turbiilcnce the events with randomly 
occurring large values of the transmission coefficient will 
be dominant in the postselected measurements. In fact, 
these are the wanted events originating from the nonclas- 
sical light source. 

Realistic radiation sources generate quantum states, 
which differ from the perfect Bell states. Important ex- 
amples are parametric down-conversion sources, which 
are often used in experimental investigations. For such 
sources we study the situation for both photon-number- 
rcsolving detectors and on/ofF detectors. The latter ap- 
pear to be very sensitive to the presence of background 
radiation and dark counts. We have analyzed this prob- 
lem for realistic parameters. Unlike in the case of photon- 
number resolving detectors, for on/ofF detectors and for 
weak turbulence the Bell parameter may become ex- 
tremely small. Increasing atmospheric turbulence may 
substantially improve this situation. 
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Appendix A: Correlation coefficient for Bell states 



Here we give the expression for the correlation coefficient E(9a,0b) for the Bell state (12), assuming that all 
detectors differ from each other. In the most general case this coefficient is given by 



E{eA,eB) (Ai) 

^ PbEs {Oa, Ob) + PH,gH, (^A, ^b) + Pv^^v^ (^A, ^b) + Phb^hb (^a, ^b) + PvB-gvB (^A, ^b) + PqEq {9 a, Ob) 
PbPb (0a,0b) +Ph^Ph^ (0a, ^b) +Vv^P^^ (0a,0b) +Phb-Ph3 (0a,0b) +Pvb^vb (0a,0b) +Pon (0a,0b) ■ 



Here p^^, p^^, p^^, p^-^, ps arc given by Eqs. (23)-(26). The coefficients Ei{9a,0b) and Pj(6'a,6'b) {i = 
{0,Ha, Va,Hb, Vbj-B}) have a different form for the case of photon-number resolving and on/off detectors. For 
photon-number resolving detectors they read as 



Eo {0A, Ob) = e-^- (iV,^ - N^^) {N^^ - N^^) , (A2) 



9a, Ob)= - + (1 - V^JN^J cos^ 9a - + (1 - VtJN^J siii^ Oa] 

X [('/t3 + (1 - '/tb)^tJ cos2 0b - (?7r3 + (1 - %b)^i^b) ^in' ^b] 
- [i^^A + (1 - VtJN^J cos2 0a - + (1 - sin' f?A] 

X [(%B + (1 - %b)^kb) cos' 0b - (%B + (1 - %b)^tb) sin' 0b] 

+ 1 [(%A + (1 - '/ta)^tJ + + (1 - %jA^aJ] 

X [(%B + (1 -^tb)^tb) + (%B + (1 -%b)^Rb)] sin 20a sin 20b cos I , 
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(A3) 



i?H,<3,(^A,0B)=e 



^^Am + (1 " 'K(B))(^T^(B) ~ ^I^A(B)) 



COS 



^A(B) 



%A(B) + (1 - '?Ha(B))(^Ha(B) - ^T^(B)) 



(A4) 



sm 



^A(B) 



i?V,,3,(^A,^B)=e-^- (iV,3<A,-^H3<,, 

X { [v^AiU, + (1 - '/t^(B))(^T^(B, - ^Ha(B,) 



(A5) 



sm 



^A(B) 



'^Hajb, + (1 - %a(b,)(^Ha(b) - ^T^(B))] cos' 0A(B)} 



Po {Oa, Ob) = e-^- (iV,^ + iV^^J (iV,^ + ^rb) > 



(A6) 



+ (1 - %a)A^ra) cos^ 0a - (77t^ + (1 - VtJN-tJ sin^ 0a] 

X [(%B + (1 - %b)^tb) cos^ 0b - (%3 + (1 - %3)iVR3) sin2 0b] 
+ (1 - %a)^ta) cos2 0a - + (1 - %a)^Ra) sin' ^a] 

X [(%3 + (1 - %B)^nB) cos' 0B - (%3 + (1 - %b)^tb) sin' 0b] 

+ ^ [(%A + (1 - '?Ta)^tJ + (%, + (1 - % J^VrJ] 

X [(%B + (1 - ??tb)A^tb) + (%B + (1 - ?7rb)^Rb)] sin 20a sin 20b cos v5 I , 



(A7) 



PH,<3,(^A,0B)-e-^- (a^T3<,,+A^R3<a,) 



^^A,U, + (1 - ''T^(3))(iVT^(3) + A^R^,3,) 



COS 



^A(B) 



A(B)^ 



(A8) 



sm 



^A(B) 



9A,^B)=e-^- (^T3,,,+A^R3,A)) 

: { [%A(B, + (1 - '7t^(B,)(^T^(B, + ^«A(B,) 



(A9) 



sm 



^A(B) 



A(B) + ^T,,3,)] COS'0A(B)} 



where Nnc is the total mean value of noise counts, cf. Eq. (38). For the case of using on/ofF detectors the coefficients 
in Eq. (|Aip are given by 



E„ (0a,0b) = e-^'- (e^-A _ e^-*) (e^^B _ e^^B) , 



(AlO) 
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= - [e^^A {r]n, + (i _ r,^J(i _ e"^^*)) cos^^a - e^-A (rj^^ + (i _ rj^J^i _ g-^-A)) sin^^A] 

X [e^-B (ry,^ + (1 - n^J{l - e-^-B )) cos^ - e^^B (rj^^ + (l - n^J{l - e'^^B )) sin^ ^b] 

- [e"^"* (%A + (1 - %a)(1 - e-"^"*)) cos2 ^a - e^-A (r;,^ + (i _ ^„ J(l _ e'^-A)) sin^ ^a] 

X [e^-B (rjn, + (1 - - e-^-B)) cos^ ^b - e^-B (^^^ + (i _ r,^J(i _ e'^-B)) sin^ ^b] 

+ 1 [e""^^ ('?T. + (1 - r?Tj(l - e-^-A)) + e^a, (^^^ + (i _ _ g-iVn,))] 

X [e^^B (^^^ + (i_^^^)(i_e-^TB))+eJVR3 (^^^ + (i _ _ e-^^B))] sin2eAsin2^BCOS(p}, 



(All) 



(A12) 



{ 



COS 



^A(B) 



sm 



^A(B) 



(A13) 



'-{ K,„e^-A(B) + (1 _ ry,^^^,)(e^-A(B, _ e^K,^,,) 



sin^ 6Ia(b) 



- [^a,(,,e^-A(B) +(i_^^^^^^)(e^H,(3) _e^-A(B,)]cos^^A(B)}, 
Po {Oa, Ob) = e-^"= (e^^* + e^^A _ 2) (e^^B + e^^B _ 2) , 



(A14) 



- e^-A (rj^^ + (i _ r,^J(i _ e'^-A)) sin^ ^a] 

(1-??Rb)( 

)) sin^^A] 



= [e^^^- (??a, + (1 - %J(1 - e-"^^-)) cos2 0A 

X [e^^B (,y^^ + (1 _ r?^^)(i _ e-^^B)) cos^ ^b - e^^e (^^^ + (i _ _ e-^^B)) sin^ ^b 

[e'^-A (^^^ + (1 _ _ e-^-A)) cos2 ^A - e^BA (r,^^ + (l _ J(l _ e 



(A15) 



COS Wb — e 



(l-r/T3)(l-e-^-B)) sin^^E 



+ 

X 



[e^-B (^^^ + (1 - ^^J(l - e ,^ 

I ('?TA + (1 - ^tJ(1 - e-^-A)) + e^H, (^^^ + (1 _ ^^J(l _ e-iVK,))] 

c [e^^B (r]^^ + (1 _ ry^^)(i _ g-^^B)) + e^^B (ry^^ + (i _ rj^J{i _ g-^^B))] sin2eAsin2^BCOS(/p}, 



(A16) 



cos^6' 



A(B) 



x{ K,3,e'^"^'^' +(l-%A<B,)(e'^"*''^' +e^-A(B, _2) 
- K^eje"^""'"' + (1 - ^RA(B))(e'^""<"' + e"^""'"' " 2)] sin^ ^a(b)} 



Pv,(e) (^A,fe)= e-N"= (e'^^BCA) + e'^'^B(A) _ 2 

+ KcBje"^"*^"' + (1 - r?K,<3,)(e'^"*™ - e"^""'"' - 2)] cos^^acb)}- 



(A17) 



x{ ^(3)6"^"*'"' + (1 - ^T,(3,)(e'^"*^"' + e"^"*^"' - 2)] sin'^ACB) 
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Appendix B: Coefflcients in the detection probabilities for the parametric down-conversion source 

Here we give the explicit form of the coefBcients Co, Ci^ , Ci^ , and Ci^^i^ in the most general case. These coefficients 
appear in the expressions for the probabilities -PiA.iB (^Aj ^b), cf. Eqs. (36) and (37), in the case of using a PDC source 

of entangled photons. They are given by 

Co = [1 + (I^Ta - 1) tanh^ x] [l + (^^Ra " 1) tanh^ x] [l + (^^Tb " 1) tanh^ x] [l + (^^Rb - 1) tanh^ x] (Bl) 

- [1 + (riTA - 1) tanh^ x] [1 + (Otb - 1) tanh^ x] ^^Ta.Tb " [l + (^^Ta - 1) tanh^ x] [l + {^k^ - 1) tanh^ x] ^Ta,Rb 

- [1 + (f^R^ - 1) tanh^ x] [1 + (^^RB - 1) tanh^ x] ^Ra,Rb - [l + (^^Ra " 1) tanh^ x] [l + (^^Tb " 1) tanh^ x] I^Ra.Tb 
+-Do, 



Cta = -f^TA [1 + (f^RA - 1) tanh^ x] [1 + (f^TB - 1) tanh^ x] [l + (^^Rb - 1) tanh^ x] (B2) 

tanh^ x [l + (^^Tb - 1) tanh^ x] £*Ta,Tb + ^^Ta tanh^ x [l + (^^Rb - 1) tanh^ x] -Dta.Rb 
+ [1 + (Qr^ - 1) tanh^ x] [1 + (Qrb - 1) tanh^ x] Z^Ra.Rb + [l + (^^Ra - 1) tanh^ x] [l + (Otb - 1) tanh^ x] Z^Ra.Tb 
--Do, 

Cr^ = -Qr^ tanh^ X [1 + (Ota - 1) tanh^ x] [l + {n^^ - 1) tanh^ x] [l + (I^Rb - 1) tanh^ x] (B3) 
+ [1 + (Qta - 1) tanh^ x] [1 + (J^Tb - 1) tanh^ x] Dta.Tb + [l + (^^Ta - 1) tanh^ x] [l + (Or^ - 1) tanh^ x] ^'ta.Rb 
+f2R^ tanh^ X [1 + (^^Rb - 1) tanh^ x] -Dr^.r^ + Qr^ tanh^ x [l + (^^Tb - 1) tanh^ x] ^'ra.Tb 
-£>o, 



Ctb = -^^Tb tanh^ X [1 + (Ota - 1) tanh^ x] [l + (Ora - 1) tanh^ x] [l + (f^RB - 1) tanh^ x] (B4) 
+f7TB tanh^ X [1 + (r!TA - 1) tanh^ x] -Dta,Tb + [l + (^^Ta " 1) tanh^ x] [l + (Orb - 1) tanh^ x] -Dta.Rb 
+ [1 + (f^R^ - 1) tanh^ x] [1 + (f^RB - 1) tanh^ x] i'RA.RB + ^^Tb tanh^ x [l + (Or^ - 1) tanh^ x] i^RA.TB 
--Do, 

Crb = -Or^ tanh^ X [1 + (^^Ta - 1) tanh^ x] [l + (^^Ra - 1) tanh^ x] [l + {^t^ - 1) tanh^ x] (B5) 

+ [1 + (f^TA - 1) tanh' x] [1 + (f^TB - 1) tanh^ x] -Dta.Tb + ^^Rb tanh' x [l + (^^Ta - 1) tanh" x] -Dta.Rb 
+r!R3 tanh^ X [1 + (^^Ra - 1) tanh^ x] -Dra,Rb + [l + (^^Ra - 1) tanh^ x] [l + (^^Tb - 1) tanh^ x] ^^Ra.Tb 
--Do, 



Cta,Tb = f^TA^^TB tanh^ X [1 + (^^Ra - 1) tanh^ x] [l + (^^RB - 1) tanh^ x] (B6) 
-Ota^^Tb tanh^ X-Dta,Tb - ^^Ta tanh^ X [l + (^Rb - 1) tanh^ x] -Dta.Rb 

- [1 + (Qr^ - 1) tanh^ x] [1 + (Or3 - 1) tanh^ x] -Dr^.r^ - Qtb tanh^ x [l + (Or^ - 1) tanh^ x] ^'ra.Tb 
+-Do, 

Cta,Rb = f^TA^^RB tanh* X [1 + (f^RA - 1) tanh^ x] [l + (f^TB - 1) tanh^ x] (B7) 
-Qta tanh^ X [l + (^^Tb - 1) tanh^ x] -Dta,Tb - ^^Ta^^Rb tanh"* X-Dta,Rb 

-Qrb tanh^ X [1 + (rjRA - 1) tanh^ x] -Dr^.Rb - [l + {^K^ - 1) tanh^ x] [l + (I^Tb " 1) tanh^ x] -Dr^,Tb 
+-Do, 

Cr^,Tb = ^Ra^Tb tanh^ X [1 + (^^Ta " 1) tanh' x] [l + (^^Rb - 1) tanh^ x] (B8) 
-riTB tanh^ X [1 + (Ota - 1) tanh^ x] -Dta.Tb " [l + (^^Ta ^ 1) tanh^ x] [l + (f^RB - 1) tanh^ x] -Dta,Rb 
-17ra tanh^ X [l + (^^Rb - 1) tanh^ x] -Dra.Rb - ^Ra^Tb tanh"' X-Dra,Tb 
+-Do, 
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Cr^.Rb = ^Ra^R^ tanh^ X [1 + {^Ta - 1) tanh^ x] [l + (Ot^ " 1) tanh^ x] (B9) 
- [1 + {Qta - 1) tanh^ x] [l + (^^Tb - 1) tanh^ x] ^Ta.Tb - ^^Rb tanh^ x [l + (^^Ta - 1) tanh^ x] ^'ta.Rb 
-^^Ra^^Rb tanh"^ X-Dra,Rb - ^Ra tanh^ X [l + (^Tb - 1) tanh^ x] -Dra.Tb 
+Do, 



where 



-Dta.Tb = ^Ra^Rb tanh^ X |e sin cos + cos sin I ^ , (BIO) 



-Dta ,Rb = ^Ra ^Tb tanh^ x | c *''' sin 6a sin 6b — cos 6a cos ^b | , (B 1 1 ) 

-Dra,Rb = ^^Ta^Tb tanh^x |~e~*'''cos^Asin^B — sin^ACOS^e]^ , (B12) 

-Dra,Tb = ^^Ta^Rb tanh^x |~e~*'''cos^ACOS0B + sin^Asin^Bj^ , (B13) 

Do = ^^Ta^^Ra^^Tb^Rb tanh^x, (BM) 

^^Ta(b, ='yTA(B,|T^A(B)|', (B15) 

^^Ra(b) =^Ra(b) I^A(B)|'- (B16) 
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